Asymptotic exchange coupling of quasi-one-dimensional excitons in carbon nanotubes 
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An analytical expression is obtained for the biexciton binding energy as a function of the inter- 
exciton distance and binding energy of constituent quasi-one-dimensional excitons in carbon nano- 
tubes. This allows one to trace biexciton energy variation and relevant non-linear absorption under 
external conditions whereby the exciton binding energy varies. The non-linear absorption lineshapes 
calculated exhibit characteristic asymmetric (Rabi) splitting as the exciton energy is tuned to the 
nearest interband plasmon resonance. These results are useful for tunable optoelectronic device 
applications of optically excited semiconducting carbon nanotubes, including the strong excitation 
regime with optical non-linearities. 
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Single-walled carbon nanotubes (CNs) — graphene 
sheets rolled-up into cylinders of ~ 1—10 nm in diameter 
and ~ lyum up to ~ 1 cm in length [l],[2j — are shown to be 
very useful as miniaturized electromechanical and chem- 
ical devices scanning probe devices [U, and nanoma- 
terials for macroscopic composites The area of their 
potential applications was recently expanded to nanopho- 
tonics [aH] after the demonstration of controllable single- 
atom incapsulation into single- walled CNs [8|], and even 
to quantum cryptography since the experimental evi- 
dence was reported for quantum correlations in the pho- 
toluminescence spectra of individual nanotubes [9( . 

The true potential of CN-based optoelectronic device 
applications lies in the ability to tune their properties in a 
precisely controllable way. In particular, optical proper- 
ties of semiconducting CNs originate from excitons, and 
may be tuned by either electrostatic doping |10| . or via 
the quantum confined Stark effect (QCSE) by means of 
an electrostatic field applied perpendicular to the CN 
axis [ll[. In both cases the exciton properties are me- 
diated by collective plasmon excitations in CNs [l2j ■ In 
the case of the perpendicularly applied electrostatic field, 
in particular, we have shown recently [ll| that the QCSE 
allows one to control the exciton-interband-plasmon cou- 
pling in individual undoped CNs and their (linear) optical 
absorption properties, accordingly. 

Here, I extend our studies to the strong (non-linear) 
excitation regime whereby photogenerated biexcitonic 
states may be formed in CNs (observed recently in single- 
walled CNs by the femtosecond transient absorption 
spectroscopy technique [HI). An analytical (universal) 
expression is obtained for the biexciton binding energy 
as a function of the inter-exciton distance and the bind- 
ing energy of constituent excitons. The formula is consis- 
tent with the numerical results reported earlier flil [l5| , 
and is advantageous in that it allows one to trace biexci- 
ton energy variation and relevant non-linear absorption, 
accordingly, as the exciton energy is tuned to the near- 



est interband plasmon resonance by means of the QCSE. 
The non-linear absorption lineshapes are calculated close 
to the first interband plasmon resonance for the semicon- 
ducting (11,0) CN (chosen as an example) under resonant 
pumping conditions [16]. They exhibit the characteristic 
asymmetric splitting behavior similar to that reported for 
the linear absorption regime [ll| . This effect could help 
identify the presence and study the properties of biexci- 
tons in individual single- walled CNs, which is not an easy 
task under non-linear excitation because of the strong 
competing exciton-exciton annihilation process [l7l - fl9l |. 

The binding energy of the biexciton in a small-diameter 
(~ 1 nm) CN can be evaluated by the method pioneered 
by Landau l20ll. Gor'kov and Pitaevski [2l[, Holstein 
and Herring |22| - - from the analysis of the asymp- 
totic exchange coupling by perturbation on the config- 
uration space wave function of the two ground-state one- 
dimensional (ID) excitons. Using the cylindrical coor- 
dinate system with the z-axis along the CN axis and 
separating out circumferential and longitudinal degrees 
of freedom of each of the excitons by transforming their 
longitudinal motion into their respective center-of-mass 
coordinates [Tl[ [23[ , one arrives at the biexciton Hamil- 
tonian of the form [see Fig. Q] (a)] 
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Here, 21,2 = 2 e i,2 — Zhi. 2 are the relative electron-hole 
motion coordinates of two ID excitons separated by the 
center-of-mass-to-center-of-mass distance AZ = Z 2 — Z\ , 
20 is the cut-off parameter of the effective (cusp-type) 
longitudinal electron-hole Coulomb potential. Equal elec- 
tron and hole effective masses m ei h are assumed [2J| and 
"atomic units" are used [2QM2^ j . whereby distance and 
energy are measured in units of the exciton Bohr radius 
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FIG. 1: (Color online) (a) Schematic of the exchange coupling 
of two ground-state ID excitons to form a biexcitonic state 
(arb. units). Two collinear axes, zi and Z2, represent indepen- 
dent relative electron-hole motions in the 1st and 2nd exciton, 
have their origins shifted by AZ, the inter-exciton center-of- 
mass separation, (b) The coupling occurs in the configuration 
space of the two independent longitudinal relative electron- 
hole motion coordinates, z\ and Z2, of each of the excitons, 
due to the tunneling of the system through the potential bar- 
riers formed by the two single-exciton cusp-type potentials 
[bottom, also in (a)], between equivalent states represented 
by the isolated two-exciton wave functions shown on the top. 



and Rydberg energy, a* B and Ry* = h 2 /(2[ia*g ), respec- 
tively, /i (~m e /2) is the exciton reduced mass. First two 
lines in Eq. (JlJ represent two non-interacting ID excitons 
with their individual potentials symmetrized to account 
for the presence of the neighbor a distance AZ away, 
as seen from the Z\- and ^-coordinate systems treated 
independently [Fig. [T] (a)]. Last two lines are the inter- 
exciton exchange Coulomb interactions — electron-hole 
and hole-hole + electron-electron, respectively [25] . Biex- 
citon binding energy is Exx = E g — 2Ex , where E g is the 
lowest eigenvalue of Eq. (TT|), Ex — —Ry*/i% is the single 
exciton binding energy with vq being the lowest-bound- 
state quantum number of the ID exciton 23[. Negative 
Exx indicates that the biexciton is stable with respect 
to dissociation into two isolated excitons. 

The Hamiltonian ([T]) is effectively two dimensional in 
the configuration space of the two independent relative 
motion coordinates, Z\ and z 2 . Figure Q] (b) , bottom, 
shows schematically the potential energy surface of the 
two closely spaced non-interacting ID excitons [second 



line of Eq. (JTJ] in the (zi, z 2 ) space. The surface has four 
symmetrical minima [representing isolated two-exciton 
states shown in Fig.JJJ(b), top], separated by the poten- 
tial barriers responsible for the tunnel exchange coupling 
between the two-exciton states in the configuration space. 
The coordinate transformation x — (z\ — Z2 — AZ) 
y = (zi + z 2 )/V2 places the origin of the new coordi- 
nate system into the intersection of the two tunnel chan- 
nels between the respective potential minima [Fig [T](b)1, 
whereby the exchange splitting formula of Refs. [2(J-[22| 
takes the form 



E g , u (AZ)~2E x = TJ{AZ), 



(2) 



where E gu are the ground-state and excited-state ener- 
gies [eigenvalues of Eq. (JTJ)] of the two coupled excitons as 
functions of their center-of-mass-to-center-of-mass sepa- 
ration, and 



9 nAZ/y/2 

J(AZ) = -J_ jhj 



-AZ/V2 



tp(x,y) 



dip(x, y) 



dx 



(3) 



x=0 



is the tunnel exchange coupling integral, where tp(x,y) is 
the solution to the Schrodingcr equation with the Hamil- 
tonian (JTJ transformed to the (x, y) coordinates. The fac- 
tor 2/3! comes from the fact that there are two equivalent 
tunnel channels in the problem, mixing three equivalent 
indistinguishable two-exciton states in the configuration 
space [one state is given by the two minima on the x-axis, 
and two more are represented by each of the minima on 
the y-axis — compare Fig.Q](a) and (b)]. 

The function ip(x,y) in Eq. ([3]) is sought in the form 

ip(x,y) = ip (x,y)exp[-S(x,y)] , (4) 

where ip = v^ 1 exp[-(|zi(x, y, AZ)\ + \z 2 (x, y, AZ)\)/v } 
is the product of two single-exciton wave functions [26( 
representing the isolated two-exciton state centered at 
the minimum z\ = z 2 = (or x = —AZ/y2, y = 0) of the 
configuration space potential [Fig. Q] (b)], and S(x,y) is 
a slowly varying function to take into account the devia- 
tion of tp from ipo due to the tunnel exchange coupling to 
another equivalent isolated two-exciton state centered at 
21= AZ, z 2 = -AZ (or x = AZ/V2, y = 0). Substituting 
Eq. ((4]) into the Schrodinger equation with the Hamilto- 
nian ([1} pre-transformed to the (x, y) coordinates, one 
obtains in the region of interest (zq dropped [2(|) 
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up to negligible terms of the order of the inter-exciton van 
der Waals energy and up to second derivatives of S. This 
equation is to be solved with the boundary condition 
S(—AZ/y/2, y) = originating from the natural require- 
ment ijj(-AZ/V2,y)=tp (-AZ/y/2,y), to result in 
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After plugging Eqs. (J5)) and Q into Eq. Q, and retain- 
ing only the leading term of the integral series expansion 
in powers of Vq subject to AZ > 1, one obtains 



J(AZ) 
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The ground state energy E g of two coupled ID excitons 
in Eq. ^ is now seen to go through the negative mini- 
mum (biexcitonic state) as AZ increases. The minimum 
occurs at AZq = isq/2, whereby the biexciton binding 
energy is E xx = - J{AZ ) = -(l/9^ 2 )(e/3) 21 '^ 1 . In ab- 
solute units, expressing vq in terms of E x , one has 



E 



xx 
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The energy Exx can be affected by the QCSE as \Ex\ 
decreases quadratically with the perpendicular electro- 
static field applied pjj. The field dependence in Eq. ([7j 
mainly comes from the pre-exponential factor. So, |-Exx| 
will be decreasing quadratically with the field, as well, for 
not too strong perpendicular fields. At the same time, 
the equilibrium inter-exciton separation in the biexciton, 
AZq = vq/2 =l/(2y / \Ex\ ) (atomic units), will be slowly 
increasing with the field, consistently with the lowering of 
I .Exx |- in the zero field, assuming \E X \ ~ r~ ' 6 (r is the 
dimensionless CN radius) as reported earlier from vari- 
ational calculations [27], one has \Exx \ ~ r~ o e as well, 
which is weaker than the r -1 dependence of Ref. 14j, 



but agrees qualitatively with the recent advanced Monte- 
Carlo simulations [HI . Interestingly, as r goes down, the 
ratio \Exx / Ex\ in Eq. (JTJ) slowly grows up approach- 
ing the ID limit l/3e « 0.12. This tendency can also be 
traced in the Monte- Carlo data of Ref. [Hj]. Finally, AZo 
goes down with decreasing r, thus explaining experimen- 
tal evidence for enhanced exciton-exciton annihilation in 
small diameter CNs [l^-Gij]- 

FigureCHshows the difference E g (AZ)-2E x =-J(AZ) 
calculated from Eqs. ([2]) and ([5]) for a specific example of 
the coupled pair of the first bright excitons in the semi- 
conducting (11,0) CN exposed to different perpendicular 
electrostatic fields. The inset shows the field dependences 
of E X x [as given by Eq. Q] and of AZ . All the curves 
are calculated using Ry* = 4.02 eV, \E X \ = 0.76 eV, and 
the field dependence of Ex reported earlier in Ref. [HI ■ 
They exhibit typical behaviors discussed above. 

Now consider the exciton absorption lineshape under 
controlled (e.g., by the QCSE) variable exciton- interband 
-plasmon coupling |ll| . In the linear (longitudinal) exci- 
tation regime, one has for the exciton with the energy e 
close to a plasmon resonance the lineshape of the form 
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J (e) [{x-ef + Aa 



[(x- 



X 2 /4] 



e) 2 (Ax 2 p + Ae 2 ) ' 



(8) 

where io=r( e )/27r, r is the spontaneous decay rate into 
plasmons, X — *J4k Ax p Iq , Ax p is the half-width-at- 
half-maximum of the plasmon resonance with the en- 
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FIG. 2: (Color online) Difference E g (AZ) - 2Ex for the cou- 
pled pair of the first bright excitons in the (11,0) CN as a 
function of the center-of-mass-to-center-of-mass inter-exciton 
separation AZ and perpendicular electrostatic field applied. 
Inset: biexciton binding energy Exx and equilibrium inter- 
exciton separation AZq (y- and i-coordinates, respectively, 
of the minima in the main figure) as functions of the field. 



ergy x p , and As is an additional exciton energy broad- 
ening (normally attributed to the exciton-phonon scat- 
tering with the relaxation time T p h). All quantities in 
Eq. are dimensionless, i.e. normalized to 270, where 
70 = 2.7 eV is the C-C overlap integral, and the condition 
assumed to hold. 
The non-linear optical susceptibility is proportional 
to the linear optical response function under resonant 
pumping conditions [lfl- This allows one to use Eq. © 
to study the non-linear excitation regime with the pho- 
toinduccd biexciton formation as the exciton energy is 
tuned to the nearest interband plasmon resonance. The 
third-order longitudinal CN susceptibility is then of the 
form [TJ, [3 



X {3) (x)= X o I(x) 
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where e xx = Exx /270 is the dimensionless binding en- 
ergy of the biexciton composed of two (ground-internal- 
state) excitons, and xo is the frequency-independent con- 
stant. The first and second terms in the brackets repre- 
sent bleaching due to the depopulation of the ground 
state and photoinduced absorption due to exciton-to- 
biexciton transitions, respectively. 

Figure [3] compares the linear response lineshape ([8]) 
with the imaginary part of Eq. ([§]) representing the non- 
linear optical response function under resonant pump- 
ing, both calculated for the first bright (ground-internal- 
state) exciton in the (11,0) CN, as its energy is tuned to 
the nearest interband plasmon resonance (vertical dashed 
line in Fig. [3]) [H|. In this calculation, Exx was taken 
to be —0.059 eV as given by Eq. ([7]) in the zero field [29| . 
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FIG. 3: (Color online) [(a), (b), and (c)] Linear (top) and non- 
linear (bottom) response functions as given by Eq. (O and by 
the imaginary part of Eq. ((9| , respectively, for the first bright 
exciton in the (11,0) CN as the exciton energy is tuned to the 
nearest interband plasmon resonance (vertical dashed line). 
Vertical lines marked as X and XX show the exciton energy 
and biexciton binding energy, respectively. Dimensionless en- 
ergy is defined as [Energy] /2'jo, where 70 = 2.7 eV is the C-C 
overlap integral. 



(Weak field dependence of Exx does not play an essential 
role here as \Exx\ "C \Ex\ = 0.76 eV regardless of the 
field strength.) The phonon relaxation time r p h = 30 fs 
was used as reported in Ref . [28[ , since this is the shortest 
one out of possible exciton relaxation processes, includ- 
ing exciton-exciton annihilation (r ee ~ 1 ps ). Rabi 
splitting ~0.1 eV is seen both in the linear and in non- 
linear excitation regime, indicating the strong exciton- 
plasmon coupling both in the single-exciton and in biex- 
citon states, almost unaffected by the phonon relaxation. 

This effect can be used in tunable optoelectronic de- 
vice applications of small-diameter semiconducting CNs 
in areas such as nanophotonics, nanoplasmonics, and cav- 
ity quantum electrodynamics, including the strong exci- 
tation regime with optical non-linearities. In the latter 
case, the experimental observation of the non-linear ab- 
sorption line splitting predicted here would help iden- 
tify the presence and study the properties of biexcitonic 
states (including biexcitons formed by excitons of differ- 
ent subbands [l3|) in individual single- walled CNs, due 
to the fact that when tuned close to a plasmon resonance 
the exciton relaxes into plasmons at a rate much greater 
than ( ^t" 1 ), totally ruling out the role of the com- 
peting exciton-exciton annihilation process. 
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